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A rate equation approach to model the denaturation or replication 


behavior of the SARS coronavirus 


J. Liu 


Abstract As a newly emerging virus, little is known about 
the SARS coronavirus, whose outbreak has brought away 
several hundred people’s lives over the world in the year of 
2003 and is seriously imperiling the human health. 
Revealing the denaturation and replication mechanisms of 
SARS coronavirus has great importance for successfully 
fighting SARS. However, experiments related to SARS 
coronavirus are extremely dangerous and therefore 
restricted only to certain specific labs with high safety 
standard. Clearly, predicting the behaviors of SARS coro- 
navirus in a wide variety of environmental conditions, 
which are not easily accessible, are thus critically neces- 
sary. In this study, we proposed to quantify the survival 
time of SARS coronavirus either in vitro or in vivo, through 
introducing thermal rate process models established from 
the well-known Arrhenius law. The complex physical and 
chemical behaviors of the SARS coronavirus can then be 
attributed to its activation energy, frequency factor, 
damage function as well as the surrounding environmental 
conditions. Based on the first data on stability and resis- 
tance of SARS coronavirus measured by members of WHO 
laboratory network, the rate coefficients involved in the 
above equations were estimated for the first time. Predic- 
tions on the survival time of SARS coronavirus in different 
temperature scale were then performed. It was found 
theoretically that, such survival time falls in an extremely 
wide range, say from several seconds in high temperature 
to an almost infinitely long time in a low temperature 
environment, which has already or is being supported by 
the currently available tests data. Applications of the 
present theory to interpret several existing phenomena 
were presented and their implementations in developing 
new technical ways for SARS prevention and clinical 
therapy were discussed. Uncertainties involved in the 
theoretical models were also analyzed and predicted. 
Parametric studies were performed to test the effects of the 
rate coefficients to the survival time of SARS coronavirus. 
Some important factors, which can significantly vary the 
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denaturation or replication process of SARS coronavirus 
were pointed out. Through regulating the parameters 
involved in the equation, certain potential therapies either 
through drug delivery or engineering approach to treat the 
SARS disease can possibly be established. Extension of the 
present model for further studies was also suggested. This 
study opens a new theoretical way for probing into the 
complex behaviors of SARS coronavirus. 


Modellierung der Denaturierung oder Replizierung 

von SARS-Korona-Viren 
Zusammenfassung Der Kenntnisstand tiber die Eigen- 
schaften des in 2003 neu aufgetretenen SARS Korona 
Virus, der einige Hundert Menschenleben gekostet hat, ist 
relativ gering. Die Ermittlung des Denaturierungs- und 
Replizierungsmechanismuses des SARS Virus ist fiir seine 
Bekaémpfung von hoher Bedeutung. Experimentelle 
Untersuchungen an diesem extrem gefahrlichen Virus 
diirfen nur durch Laboratorien mit einem hohen Sicher- 
heitsstandard erfolgen. Die Vorhersage des Verhaltens des 
SARS Virus in unterschiedlichen Umgebungsbedingungen 
ist dabei erforderlich. In der vorliegenden Studie wird die 
Uberlebensdauer des Virus unter Labor- und realen 
Bedingungen durch Anwendung der bekannten Arrhe- 
nius-Beziehung fiir temperaturabhangige Vorgange 
ermittelt. Das physikalische und chemische Verhalten des 
SARS Virus wird anhand der zugrundeliegenden Modell- 
Parameter beschrieben. Basierend auf den ersten Mes- 
sungen von Mitgliedern des WHO-laboratory-network 
tiber die Stabilitat und Widerstandsfahigkeit des Virus 
wurden erstmalig die Geschwindigkeitskoeffizienten des 
Berechnungsmodells bestimmt. Vorhersagen der Uberle- 
bensdauer des SARS-Virus unter unterschiedlichen Tem- 
peraturbedingungen wurden ausgefihrt. Das sich hieraus 
ergebende, sehr unterschiedliche Ausmaf der Uberle- 
bensfahigkeit in Abhangigkeit der Umgebungstemperatur 
ist durch den Vergleich mit verfiigbaren experimentellen 
Ergebnissen bestatigt worden. Die Anwendung der vor- 
gestellten Modellierung zur Interpretation realer Phano- 
mene und zur Entwicklung technischer Mafnahmen zur 
Vorbeugung und klinischen Therapierung von SARS wird 
diskutiert. Der Einfluf von Unsicherheiten des Modells 
wird analysiert und abgeschatzt. Parametrische Studien 
sind durchgefiihrt worden, um den Einfluf der Ge- 
schwindigkeitskoeffizienten auf die Uberlebensdauer des 
SARS Virus darzustellen. Einige wichtige EinfluSgréfen 
auf die Denaturierung und Replikationsfahigkeit des SARS 
Virus werden aufgezeigt. Durch eine Variation der Mo- 
dellparameter kann die potentielle Wirksamkeit medika- 
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mentéser oder physikalischer Therapien abgeschatzt wer- 
den. Erweiterungsméglichkeiten des vorgestellten Modells 
werden vorgeschlagen. Die vorliegende Studie erméglicht 
neue, theoretische Vorgehensweisen zur Untersuchung des 
komplexen Verhaltensmusters des SARS Virus. 


List of symbols 

A surface area of droplet 

C specific heat of droplet 

C(0) virus concentrations at initial state 
C(t) virus concentrations at time t 


Eq activation energy (J/mol) 

Eac activation energy of virus 

[Eo] total concentration of enzyme 

[E,] | concentration of active enzyme 

[Ei] concentration of denatured enzyme 
Fy injury fraction 

h heat convection coefficient 

K chemical reaction rate 

K, rate constant for immune reaction 
K, reaction equilibrium constant 

Kg equilibrium constant 

P collision number or frequency factor (s~') 
P, frequency factor of virus (s~') 

R constant of perfect gas (8.31 J/mol/K) 
[S] virus concentrations 

[S*] concentration of all protein substances 
[So] protein concentration 

t time (s) 

to survival time (s) 

T absolute temperature (K) 

Vv volume of droplet 

Vi reaction rate 

V2 reaction rate 


Greek symbols 

ratio for volume of all protein substances to 
that of SARS coronavirus 

p density of droplet 

Q injury function 

A uncertainty of variable 

AH _ variation of enthalpy 
A 


R 


S variation of entropy 
Subscripts 
crit critical value 
i initial state 
opt optimum 
oe) surrounding air 
1 
Introduction 


Named by the world health organization (WHO), severe 
acute respiratory syndrome-SARS is a newly emerging 
world-wide infectious disease [1-3], which is also com- 
monly called atypical pneumonia in China [1-4]. Since its 
outbreak, SARS has brought about serious threatening on 
the human health. Up to now, investigating the patho- 
genesis and treatment therapy on SARS keeps receiving 


great attentions over the world. Some encouraging pro- 
gresses have been achieved in the SARS research such as 
identification of cause of disease [5-7], gene sequencing 
[8-10], virus detection [11], epidemiology analysis [12, 
13], drug and bacterin development [14], pathology [15], 
and clinical therapy [16, 17] etc. It now becomes clear that, 
the pathogeny to induce SARS is a brand new type of 
coronavirus, whose behavior turns out rather unique and 
had never been found before in nature. However, since 
little is known about the properties of SARS and no 
effective drugs are available, successful treatments on this 
disease at the present stage are extremely difficult, if not 
impossible. Clearly, revealing the physical and chemical 
behaviors of the SARS coronavirus either in or outside the 
human body will be very beneficial for fighting the SARS. 

One of the latest advancements made by the WHO 
multicentre collaborative network (MCN) on SARS is the 
successful test on the survival time of SARS coronavirus 
under several different environmental conditions [18]. 
Owing to this contribution, we now get to know that the 
survival time for the SARS coronavirus in the air spans 
from hours to days or even longer time. Such a huge dif- 
ference in time scale is rather perplexing, whose inter- 
pretation is waiting for great efforts. In fact, the survival 
behavior of SARS coronavirus is a complex issue, which 
depends both on the surroundings and the virus itself [19]. 
To reveal the intrinsic mechanisms, a systematic theory 
concise in concept would be very desirable. In this study, 
we proposed to quantify the denaturation or replication 
process of the SARS coronavirus through a rate equation 
theory and new models capable of comprehensively 
characterizing the physical and chemical behaviors of 
SARS coronavirus will be established. 


2 

Rate process model for survival time of SARS coronavirus 
It is well known that, virus is a non-cellular organism, who 
realizes its surviving by attacking a hosting cell. The virus 
could not divide itself automatically [20]. As an organism 
consisted by certain special nucleic acid and protein, the 
denaturation of SARS coronavirus can be regarded as a 
temperature dependent rate problem, just as any natural 
substances do. In fact, it has long been a basic way to 
address the thermal injury of biological tissues by using 
the rate equation. Henriques and Moritz [21] were the first 
to characterize the thermal injury as a first order chemical 
reaction process. The concise equation given by them late 
became the start point to analyze the thermal injury pro- 
blem of biological tissues. Xu and Qian [22] made new 
progress in better understanding the thermal injury 
mechanisms by introducing the denaturation of enzyme 
protein. Recently, Bhowmick and Bischof [23] performed 
an experimental test on the injury problem of the cell 
membrane. Typical biophysical and biochemical events 
occurring in thermal injury at macro or cellular level 
generally include [24]: thermal denaturation of proteins, 
alterations in metabolic process, thermally induced alter- 
ations in physical or chemical characteristics of cells like 
hyperpermeability of membrane, intracellular ionic con- 
centration and nuclear degradation, loss of birefringence 
properties in muscle and collagen, loss of hemoglobin 
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from red blood cells. Similarly, the denaturation process of 
the SARS coronavirus in the air or its hosting cell should 
also follow the same mechanism as above, although may 
vary in the rate coefficients. The rate process model is 
originally established from the chemical reaction dynamics 
[25-28]. Referring to Fig. 1, one can see that for a typical 
reaction, the thermally activated reactant should overcome 
an energy barrier E, in order to generate products. Cor- 
respondingly, possibilities for the virus to replicate or die 
out also depends on the magnitude of its stored energy, 
and the premise for this to happen is that the stored 
energy should be higher than that of its critical reaction 
value, i.e. the activation energy. If treating the denatur- 
ation process of the virus as a first order process and 
defining the injury function as ©, one can write its 
changing rate as, 


Q = dQ/dt. (1) 


Integration of the above equation over the whole period 
of thermal denaturation is Q), i.e. 


Q= [taasanae (2) 


Then, variation of the virus appears as a temperature 
dependent rate process whose chemical reaction rate K 
(i.e. injury rate) can be described by the well-known 
Arrhenius equation [25]: 


dQ K=P E, ‘ 
dt EO RT (3) 
where, R is constant of the perfect gas (8.31 J/mol/K), Eq 
the activation energy (J/mol),P the collision number or 
commonly called as frequency factor (s~'), T the absolute 
temperature (K), and ¢ the time; Q stands for injury ratio, 
when Q = 1, it is generally regarded as that the object has 
almost entirely denatured (or injured) [27, 28]. Generally, 
for different biological organisms or environmental con- 
ditions, the magnitudes for P and E, could vary in a very 
large range. 

Combining (2) and (3), one can obtain the Henriques 
equation to quantify the thermal injury of the biological 
material [21], i.e. 


= f - exp(—E,/RT)dt. (4) 


In this equation, only two coefficients E, and P are 
needed to perform a comprehensive analysis on the de- 
naturation problem of the biological samples. Such inte- 
gration equation is very easy to quantitatively characterize 
the denaturation of virus. 


reactants 


\ et 


products 


energy 


Fig. 1. Energy state diagram of a biochemical process 


Here, the physical meaning for Q can be expressed as 
the logarithm of the ratio of the concentration of the initial 
virus particles and that of the surviving virus particles, i.e. 


Q(t) = In (a3) (5) 


where, C ((0) and C (t) are virus concentrations at initial 
state and that at time t¢, respectively. 

Its meaning can further be interpreted as follows. 
Assuming the thermal process as a first order activity, one 
has [26] 
dC(t) 

—— = —-KC(t). 6 
1 (t) (6) 

For time dependent rate parameter K (due to variation 

of temperature), the virus concentration should read as 


C(t) = Coe io KO", (7) 


Therefore, it can be seen that the thermal injury in fact 
accords with (4), i.e. 


Q(t) =1n(Se) = [ xae= [ ?-20(-B)at 
(8) 


Then the injury fraction can be described as [26] 


Fy= ae =1-¢™". (9) 


The above dynamics model indicates that, if the tem- 
perature is lower than certain critical value, the rate of the 
injury accumulation could be omitted. On the contrary, 
the injury rate will significantly increase. This critical 
temperature is generally defined as that when dQ/dt 
becomes equal to 1.0, i.e. [25] 


E, 
Terit Sy ae 
Rln(P) 
(10) 
In other words, if the temperature is higher than this 


value, the thermal environment can easily kill the SARS 
coronavirus. 


Q 
a Pexp(—E,/RTerit) =1 ie. 


3 

Equation for in vitro survival time of SARS coronavirus 
The above equation is very useful in characterizing the 
survival time of SARS coronavirus. As will be shown later, 
it can provide predictions close to the reality and help get a 
clear and concise view on the complex behaviors of SARS 
coronavirus. With an extremely small size (in micro/nano 
scale), the SARS coronavirus would quickly reach a ther- 
mal equalibrium with the environment, which can be 
illustrated as follows. Assuming that a liquid droplet (such 
as sputum) containing virus was dropped to the sur- 
rounding air with temperature of T.., its temperature 
response can be described by a lumped system model, i.e. 


—hA(T — T.) = pVC— 11 
(T-T..) = pes (11) 
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where, p, C are respectively the density and specific heat of 
the droplet, A, V the surface area and volume of the dro- 
plet, h the heat convection coefficient between the droplet 
and the surroundings. 

Solution to (11) can easily be obtained as: 


T = T,. + [T; — Too] exp[—(hA/pVC)t] 


where, T; is initial temperature of the droplet. 

For an extremely small liquid droplet, the second term at 
the right hand side of (12) will become zero within seconds. 
Compared with the long survival time of the SARS corona- 
virus, which is generally in minutes or even days, the 
temperature of the SARS coronavirus can be treated as a 
constant without causing evident error. Therefore integra- 
tion of (4) will lead to the following theoretical equation for 
the survival time tf) of SARS coronavirus in the air as, 


ty = Q- exp(E,/RT)/P. (13) 


It indicates that, the SARS coronavirus is very sensitive to 
the temperature such that the survival time fp exponentially 
depends on the temperature T, which accords well with 
the currently available observation. Eq. (13) still reveals that, 
the higher activation energy or the lower temperature or 
frequency factor for the virus, the longer time for it to live. 

Should be pointed out that, application of (13) often 
implies approximations due to uncertainty of the rate 
coefficients. The reason mainly comes from the experi- 
ments, which may have a low testing accuracy or just 
record insufficient information, especially because defini- 
tion of some important values such as critical survival rate, 
and details on the substrate are only approximately true. 
Besides, experimental tests up to several days could not 
guarantee a constant environmental temperature. These all 
contribute to the statistical errors and may lead to 
inaccurate estimation on the rate coefficients. To test the 
effects from each of the parameters to the survival time in 
(13), to was expressed as the following function 


to = f (Ea; P, T,Q). (14) 


Assuming At is uncertainty caused by the approximate 
parameters, then one has 


At 


(12) 


_ Oto * 1 Oto n (Oto 2 (Oto ? 
: i (<2 ap,) +(Mtar) +(Har) + (28.0) 
(15) 


where, A stands for uncertainty for each of the variables. 
Eq. (15) can further be expressed as [29]: 


Therefore, if uncertainties for the four parameters E,, P, T 
and Q are given, their effects on the survival time can be 
evaluated and its uncertainty limit will thus be predicted. 


4 
Equation for in vivo survival behavior of SARS coronavirus 
The discussion given above is only useful when the 

SARS coronavirus lives in the air. For the process that the 
virus attacks its hosting cell, then replicates or dies, more 
factors should be taken into concern. Analysis performed 
below will be very beneficial to better understand the ef- 
fects of the concentrations of the virus and immune cell etc 
to the infection mechanisms of SARS disease. Generally, 
when virus enters the lung, it would try to replicate itself 
by attacking the hosting cell. Without adopting effective 
measure to prevent this from happening, the virus con- 
centration would become higher and higher until cause 
injury to the lung tissues [12]. Therefore, quantitative 
evaluation on the early damage of the lung and then pre- 
dict the damage time will be very important for a suc- 
cessful treatment planning on SARS. Clearly, an 
appropriate theoretical model would be desirable for such 
purpose. In the following, a theoretical model will be 
established to characterize the replication or death 
behavior of the SARS coronavirus in human body. 

As is known, replication of virus must recur to the 
enzyme of its hosting cell. Assuming this enzyme-cata- 
lyzed process follows a first order equation, then its 
reaction rate V; is proportional to the virus concentrations 
[S] and that of the active enzyme [E,], i.e. [22] 


Vi = K{S][Ea]- (18) 


Within the scale of the physiological temperature, the 
active enzyme |E,| and the denatured enzyme [Ej] obeys 
the following equilibrium relation: 


Ka 
[Ea] — [Ei] (19) 
where the equilibrium constant Kg depends on temperature 
and can be determined via the van’t Hoff equation [22], 
Ka = [Ei]/[Ea] = exp(AS/R — AH/RT) (20) 
where AH and AS are respectively the variations of the enth- 
alpy and entropy during the enzyme denaturation process. 


Since the total amount of the active enzyme and the non- 
active enzyme is a constant, i.e. [E,] + [Ei] = [Eo], one has 


[Ea] = [Eo]/[1 + Kal. (21) 


Then, the chemical reaction rate of virus replication 
under enzyme catalysis was obtained as 


At | (@lnty AE,\* | (OlntyAP\* (dint, AT\* 
tf \V\OlnE, E,/ \OmP PP) '\OlnTT/) | 


Substituting (13) into (16) leads to 
AP\*  (AQY?* 
P)/) ‘\Q)° 


*..| (BY (2) +(4)]: 


gi 
Oln to a) . (16) 


OlnQ QO 


V, = K[Eo|(S|/[1 + exp(AS/R — AH/RT)]. (22) 


It indicates that, there exists an optimum temperature 
Topt where the reaction rate is the highest, and its value can 


4 


be determined by solving equation de = 0. Generally, 


T=Tp 
this value is close to the normal physiological temperature 
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37 C. If substituting T,,; into (22), one can further obtain 
the optimum reaction rate Vop;. However, expressions for 
both Top; and Vo, will not be given here for brief. 

The enzyme catalysis would help to realize a virus 
replication process. In contrast to this, phagocytosis by the 
immune reaction due to body response will restrain the 
growth of the virus, although the detail is not clear pre- 
sently. For a generalized purpose, assuming that the cell 
concentration to induce a first order immune reaction is 
[C], the reaction rate will then also be proportional to the 
concentrations of the virus and the immune cells, i.e. 


V2 = KA[SIIC] (23) 


here, K, is the rate constant for the immune reaction, [C] 
can be regarded as a transient value and follows the rela- 
tion as: 


K, 
[¢]| —_ Iq] 
where, the reaction equilibrium constant K, can be 
approximately described by the Arrhenius equation, 
K, = P, exp(—E,-./RT) (25) 


where, E,,, P, are respectively the activation energy and 
frequency factor of virus subject to immune reaction. 
Concentration of the immune cells can be expressed as 


d[C] 


(24) 


ao = Kd (26) 
In{C]/[Co] = — | ‘Rica eR. 
(27) 


Since mechanism for enzyme-enabled reaction 
differs from that of the phagocytosis reaction, the 
optimum temperatures for the two processes are also 
different. If certain external heating can be applied to 
enhance phagocytosis while restrain enzyme enabled 
reaction, then a useful way to kill the SARS virus can 
seems be established. 

With the above relations, one can predict the survival 
time of the SARS coronavirus in the human body. Assuming 
the initial concentration of virus invading the human body 
is [So], and the rate for the enzyme enabled reaction is Vj, the 
rate for the phagocytosis reaction induced by immune cell is 
V2, then the transient virus concentration [S] can be 
determined by the following equation, 


d[s]/dt = (Vi — V2) 
K|Eo] 
= |S K2|C] > (28 
| if + exp(AS/R — AH/RT) 2 ih ) 
where the transient value of [C] is determined from (27). 
Clearly, if temperature and the concentration of 


immune cell [C] are both constant, the virus concentration 
can be obtained as 


KE] 
TnISiAI0) = oF + exp(AS/R — AH/RT) Kala}. 
(29) 


If most of the virus denatured, i.e. the current concen- 
tration [S ] significantly deviate from its initial value [S o], 
then the damage can be regarded to happen. From (5), one 
can find that the damage function for the virus reads as 
Q = In{So]/[S]. Then the survival time for the virus to 
invade the human body is obtained as 


K|Eo| = 
aise} ew) 


If parameters involved in the above equation can be 
estimated through experiment, one can predict the 
behavior of the SARS coronavirus and thus the SARS 
pathology development in the human body. However, due 
to difficulty to carry out an in vivo test, there is currently a 
strong lack of experimental data for the survival time of 
the SARS coronavirus in human body. This study still 
could not provide a calculation on such situation. 

Similar as above, due to uncertainty of the interior 
environment of human body such as variations in pH 
value, temperature, drug concentration, body response, 
substance variation even experimental measurements, 
uncertainty should also be considered when applying the 
above equation to predict the behavior of SARS corona- 
virus to attack the hosting cells. For such in vivo situation, 
the survival time can be expressed as [29], 


to = af KIC 


to =f (Eos Bas Pi Pa TQ AS AA IG, Cy) (31) 


Combination of (30) and (31) will also lead to a uncer- 
tainty equation similar to (17). But it will not be presented 
here for brief. 

Clearly, uncertain factors for the survival time of the 
SARS coronavirus in vivo appear much complex than that 
of in the air. This makes it hard to accurately predict the 
disease development. Such situation should be paid with 
special attentions in the near future. 


5 

Estimation of rate process coefficients 

of SARS coronavirus 

Before one can perform a prediction on the survival time of 
the SARS coronavirus, the rate coefficients P and E, must be 
obtained. Up to now, only very limited parameters for the 
biological samples are reported [25, 27]. Especially, there is a 
strong lack of such data for the bio-samples in cell level. 
Even for ordinary virus, no rate coefficients are available, let 
alone the corresponding information for the brand new 
SARS coronavirus. Generally, estimation of E, and P can be 
performed by systematically measuring the survival time of 
virus subject to different temperatures. Then based on a 
large quantity of the data of f) versus T, E, and P can be 
calculated through a least square fitting on the (13). The 
procedure can be as follows: placed the virus samples in a 
specific environment, and applied it with surface heating or 
cooling and kept the virus for a specific time, then complete 
the viability test. The same experiments were repeated many 
times under different constant temperatures. The critical 
damage degree in each case was chosen from a series of 
damaged virus samples, i.e Q = 1, or C(to) = 36.8%C(0). 
Then through rewriting (13) to its logarithm form, one has 
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inte) = _ inp (32) 


From the intercept between In(to) and the critical value 
of 1/T, one can get P; while from the slope of the curve, E, 
can be obtained. 

In practice, due to high infection of SARS coronavirus, 
experiments are restricted only to certain labs with high 
safety standard. Fortunately, the WHO MCN for SARS has 
experimentally obtained the first data on stability and 
resistance of SARS coronavirus, which makes the present 
study possible. Here, we will make an analysis on their 
data as modified in Table 1, where the results were par- 
ticularly chosen from those in neutral environment. 
Although these data may contain approximation, it does 
provide great valuable knowledge on the survival behav- 
iors of the SARS coronavirus in various environments. 

Based on the data in Table 1, the curve for the 
Inf) ~ 1/T can be drawn as Fig. 2, where different symbols 
represent data from different labs. It shows that, these data 
appear as some what a little disperse, which indicates that 
the experimental test may contain uncertainty. Theoreti- 
cally (see Eq. (32)), relation between In t, and 1/T follows a 
linear equation. However, due to approximation in the 
measured data, this relation does not strictly hold. 
Therefore, the presently estimated rate coefficients are in 
fact only mean values. Through averaging the data from 
different labs, we obtain the straight line in Fig. 2, which is 
just the rate equation for the survival time of the SARS 
coronavirus in the atmospheric exvironment, i,e, 


In fo = 11199.5/T — 26.4476 (33) 


Comparing between (32) and (33), one can obtain for 
the first time the activation energy and the frequency 
factor for the SARS coronavirus respectively as follows, 


E, =11199.5 x 8.31 = 0.9307 x 10°(J/mol), 
P = 3.429 x 10'!(s~'). (34) 


The above values appear much smaller than that of the 
ordinary bio-samples. The fact for the small activation 
energy implied that it is easy for the SARS coronavirus to 
react to its contacting objects, which also explain the 
phenomena that this virus can easily be destroyed by an 
ordinary disinfector [18]. Meanwhile, it also indicates that 
the SARS coronavirus has a high infection tendency. But if 
its replication path can be successfully blocked off, the 
SARS coronavirus could also not survive for a long time. 
Effects of the frequency factor appear opposite to the 
above. 


6 
Results and discussion 


6.1 

Prediction of survival time of SARS coronavirus 

Having obtained the above rate coefficients, one can then 
perform a theoretical prediction on the survival time of the 
SARS coronavirus in a wide variety of temperature ranges. 
Depicted in Fig. 3 are the predicted survival times of the 


Table 1. Experimental conditions and survival time of SARS coronavirus (from [18]) 


Lab* Substrate Condition Survival time 
GVU Virus spiked in baby stool Ph 6 - 7 (regarded here 3h 
at room temperature) 
QMH Stool Room temperature At least 2 days 
Urine Room temperature At least 24 h 
Virus culture On plastic surface At least 2 days 
medium+1% bovine serum in Room temperature 
Virus culture 30 - 37°C At least 1 h 
medium+1% bovine serum 
Virus culture 56°C Degration of titre over time (1000 
medium+1% fetal calf serum infections virus units in 15 min) 
NIID Virus culture +2% 4°C At least 4 days 
fetal calf serum 
Virus culture +2% 37°C Less than 4 days 
fetal calf serum 
Virus culture +2% 56°C Less than 30 min 
fetal calf serum 
UnivM Virus culture 4°C At least 21 days 
CUHK “virus in phosphate Room temperature on PBS Stool 
buffered saline (PBS) Plastered wall 24h 36h 
“virus in sterilized stool Plastic surface 36h 72h 
Formica surface 36h 36h 
Stainless steel 36h 72h 
Wood 12h 24h 
Cotton cloth 12h 24h 
Pig skin > 24h > 24h 
Glass slide 72h 96h 


*CUHK: Chinese University Hong Kong 
GVU: Government Virus Unit, Dept. of Health, Hong Kong, SAR China 
QMH: Queen Mary Hospital, The University of Hong Kong, Hong Kong, SAR China 


NUD: National Institute of Infectious Diseases, Tokyo, Japan 


UnivM: University Marburg, Germany 
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Fig. 2. Tests data on absolute temperature and survival time of 
SARS coronavirus (from WHO measurements [18]) 


SARS coronavirus in the temperature of 0° C-80°C. It 
indicates that, the SARS coronavirus can survive in an 
extremely wide range of time, say, from several minutes to 
hundreds of hours. Particularly, at a low temperature, 
SARS coronavirus can even live much longer. For an 
environment at temperature close to 0°C (such as in a 
refrigerator), the survival time of the SARS coronavirus 
can be as long as 540 h. While at a lower temperature such 
as -80°C, this time will sharply increase to about 

5.375 x 10° days. In other words, the SARS coronavirus 
may live forever. Therefore, as a newly emerging virus, it is 
not realistic to expect that the SARS coronavirus would 
soon completely disappear, since the virus could easily 
remain alive through the winter and then resuscitate in a 
warmer weather via various ways such that drying the 
SARS coronavirus in summer may help it realize a safe 
dormancy and then preserve itself from being thermally 
killed. Since the survived SARS coronavirus may induce a 
new cycle of SARS epidemic, the situations to confront the 
SARS disease are still very serious unless bacterin will be 
found soon. Therefore, establishing certain effective way to 
prevent SARS in advance in the season of its high infection 
is very necessary. From Fig. 3, one can observe that, the 
SARS coronavirus is very sensitive to the temperature. For 
example, at the room temperature around 25°C which is 
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Fig. 3. Temperature dependent survival time of SARS coronavi- 
rus 


comfortable for human body, the SARS coronavirus can 
live up to 16 h. Such time scale has already been long 
enough for the virus to spread to a large area of regions. 
Clearly, temperature over 20°C has produced certain effect 
to decrease the viability of the virus. To show its detail, we 
redraw the survival time curve at temperatures above 20°C 
in Fig. 3 as Fig. 4 and define the temperature above 20°C 
as thermally important temperature. In the following, we 
will make a few discussions on the survival behaviors of 
SARS coronavirus at several special temperatures. If the 
environmental temperature is at 37°C, Fig. 4 shows that 
the survival time for the SARS coronavirus is to = 3.97 h. 
Therefore, at a hot season such as summer, it is generally 
not easy for the SARS coronavirus to survive, and possi- 
bility for inducing the disease epidemic will be signifi- 
cantly decreased. On the other hand, if intentionally 
heating the virus by using the tumor hyperthermia tem- 
perature 42°C, which is bearable by human body, the 
predicted survival time for the SARS coronavirus is 

to = 2.24 h (see Fig. 4). Consequently, treatment of SARS 
through high temperature heating is possible, as long as 
the heating can be kept to a long enough time. Because 
hyperthermia therapy is a relatively safe and non-noxious 
way, it is worthy of trying to cure SARS, especially since no 
miracle drug and definite treatment strategy are available 
at the present time. Besides, if the temperature was further 
increased to a higher however still bearable level, the 
treatment time can be evidently shortened. For example, if 
choosing heating temperature 56°C to do calculation, the 
predicted survival time for the SARS coronavirus is only 
to = 0.49 h. In addition, if certain assistant drug was 
designed to decrease the activation energy or increase the 
frequency factor of the SARS coronavirus, the virus could 
possibly be killed by an ordinary heating under a relatively 
low temperature. In this sense, combining drug delivery 
and hyperthermia to cure SARS is worth of investigation. 
We have proposed an engineering approach to restrain the 
replication of SARS coronavirus through directly venti- 
lating hot oxygen to the human alveolus [30]. Should be 
pointed out that, the above analysis is performed based on 
the data of SARS coronavirus in air. Factors in human 
body are rather complex such as that there exist both 
factors favorable for immune or drug to react and that for 
the virus to replicate. Therefore, when developing a 
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Fig. 4. Thermally important temperature scale to kill SARS cor- 
onavirus 
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treatment therapy in clinics, such comprehensive factors 
should be fully considered. In certain extreme situation 
such as in a high temperature environment, the SARS 
coronavirus can only survive for a short period. From 
Fig. 4, when the temperature approaches 80°C, the pre- 
dicted result shows that the SARS coronavirus will dena- 
ture within 3 min. This conclusion has great importance 
for practice. We have proposed to heat a polluted room to 
completely kill the virus inside [31]. It seems now clear 
that, 80°C is an appropriate temperature to efficiently kill 
the virus. Killing the SARS coronavirus to clean the room 
is a new way of disinfection and produces no pollution to 
the environment. Such method is worth of investigation. 
Overall, the above results predicted from (33) accord well 
with the currently available measurement and observed 
phenomena [18]. Although the theoretical equation is 
founded from limited data, its efficiency in prediction is 
satisfactory. Theoretically, the time for the SARS corona- 
virus to live falls in an extremely wide range, but not just 
the currently reported value. This time depends on the 
magnitudes of the P, E, and T and the physical and 
chemical factors capable of changing these parameters. 
However, we should point out that, the experimental tests 
are still far from being enough, since different environ- 
mental factors and substrates will have effects on the 
survival time of the virus. Therefore, to completely 
understand the survival behavior of the SARS coronavirus, 
strictly defining the test conditions and perform corre- 
sponding work are still very necessary. 


6.2 
Effects of the rate coefficients 
In the following, we will perform an analysis on the 
uncertainty of the survival time predicted from (33). For 
this purpose, uncertainties for the activation energy, 
temperature, frequency factor and degree of denature are 
all prescribed as 10%. Then the uncertainty of the pre- 
dicted survival time at room temperature 25°C is 
At/ty = 0.88. Presented in Fig. 5 are estimated survival 
times for more situations. From that, one can find that the 
temperature and activation energy affect most to accu- 
rately predict the survival time. From (17), we can draw a 
conclusion that at a high temperature, the accuracy for 
predicting the survival time is also high, while at low 
temperature, the prediction accuracy will become worse or 
even have several times of errors. What appear more 
reliable from the theoretical prediction is the average 
survival time which overall reflects the reality. However, 
the real value usually fluctuates around this average 
magnitude. This conclusion would be very useful in cor- 
rectly interpreting the current or future measured results. 
When in the air or in human body, if certain disinfector, 
acid and alkali substance, drug, bacterin, or some chemical 
and physical factors may induce variations of the activa- 
tion energy and frequency of the SARS coronavirus, it is 
possible a simple heating process would kill the virus. On 
the contrary, the survival and replication capability of the 
SARS coronavirus may be significantly enhanced, its 
ability of anti-drug, heat resistance and infection will 
increase. Theoretically speaking, all these situations are 
possible. Parametric studies given in Fig. 6 show that, for 
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Fig. 5. Uncertainty of predicted survival time (At/f)) under dif- 
ferent temperature level and subject to uncertain rate coefficients 


different rate coefficients, survival time for the SARS 
coronavirus can vary in orders even under the same 
environmental temperature. Therefore, if one wishes to 
comprehensively analyze or test the survival behavior of 
the SARS coronavirus, the environmental conditions 
should be strictly defined. 


6.3 

Effects of enzyme denaturation and immune reaction 
From Eq. (30), one can find that, one of the most important 
factors to enable virus replication and hyperplasia is the 
enzyme. When a small part of the enzymen becomes 
irreversibly denatured, it may make it impossible for the 
SARS coronavirus to replicate. Generally, denaturation of 
enzyme is related to the temperature, pH value and drug 
application etc. If certain approaches were taken to 
intervene these factors, new therapy on the SARS can pos- 
sibly be established. On the other hand, the phagocytosis on 
the virus due to immune regulation from the human body 
will also help to kill the virus. In fact, this is one of the most 
important factors for many SARS patients to survive. 
Presently, due to lack of specific medicine, the clinical 
treatment on SARS is still above all a kind of supporting 
therapy, whose essence is to regulate the immune reaction of 
the human body and thus finally destroy the SARS coro- 
navirus. However, if inappropriate regulation was made, it 
may even enhance the replication of the virus and then make 
situation more serious. Therefore, study on the mechanisms 
of the SARS occurrence and development should be paid 
with much more attentions. 


6.4 

Denaturation of SARS coronavirus and its affinity 

with human body 

It becomes gradually clear now that it is very easy for the 
SARS coronavirus to deposit in the areas of the respiratory 
tract, alimentary canal and neural elements, which is 
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perhaps one of the reasons for such virus to have a strong 
infection. As frequently observed, the SARS coronavirus 
tends to adhere to these tissues. Therefore, revealing the 
characteristics of the physiological solutions in these areas 
will help to better understand the attacking behaviors of 
the SARS coronavirus. Similarly, we can also explain the 
other SARS phenomena. For example, the reason leads to 
virus aberrance is perhaps because its rate process has 
been changed forever. For example, some studies show 
that, the infection ability for the filial generation becomes 
weak, which may reflect that the process has changed the 
rate coefficients such as P and E,. It looks the filial gene- 
ration does have a decreased activation energy or a pos- 
sibly increased frequency factor. However, the SARS 
coronavirus does can change to a converse aberrance. This 


vival time of SARS coronavirus 


is because an extremely small change in its structure will 
possibly vary its activation energy. When the above 
parameters appear to induce an increased viability of the 
SARS coronavirus, its survival time will be lengthened and 
the infectivity will be improved. This should be born in 
mind with caution and effective prevention measure 
should be made. 


6.5 

Effects of acid-alkali environment and sterile liquid 

The rate equation established above has only considered 
the effect of a single temperature. And the estimated rate 
coefficients are obtained by using the data measured in 
neutral environment. More complex environmental situ- 
ations were not taken into account. In the coming future, it 
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is worthwhile to establish a similar rate equation, which 
however includes more factors such as pH value etc. Then 
through a limited tests on the survival time of the SARS 
coronavirus subject to different thermal as well as pH 
conditions, the corresponding rate coefficients can be 
estimated which would be beneficial for predicting more 
wide behaviors of the SARS coronavirus. Such effort would 
be of significance for the clinical treatment as well as 
prevention of SARS. For example, if one could change the 
acid and alkali conditions at the human tissues easily 
attacked by the SARS coronavirus, it will help to enhance 
the killing performance of the delivered drugs and therapy. 
For the prevention, it would be beneficial to select certain 
disinfector, which is more environmentally friendly and 
causes less pollution. 


6.6 

Cause of prolonged survival time of SARS coronavirus 

in an agglomerate with protein substances 

The ability for the SARS coronavirus to survive in the air 
significantly depends on the surroundings. However, if 
virus was contained in a protein environment such as 
saliva or respiratory secretion, its survival time will 
become much longer [32]. As is well known, the phlegm or 
saliva consists of a lot of protein and nutrimental sub- 
stance. One of the reasons for the SARS coronavirus to live 
longer in such situations is perhaps because these com- 
pounds tend to maintain large activation energy E, or a 
low frequency factor P. As predicted in Fig. 6, if only these 
two parameters change a little, they will be significant 
enough to induce a variation in orders of the survival time 
of the SARS coronavirus. 

On the other hand, the reason for the SARS coronavirus 
to survive longer can be regarded as a kind of concen- 
tration damage. If treating the SARS coronavirus deposited 
in the phlegm or saliva as also a part of this protein like 
agglomerate, then only when all the protein substances in 
the agglomerate were denatured can the packed SARS 
coronavirus be regarded as completely destroyed. In this 
way, one can predict the prolonged survival time of the 
SARS coronavirus deposited in phlegm or saliva. Substi- 
tuting (5) into (13), one obtains, 


fo = In((S"]/(S]) - exp(E./RT)/P (35) 


Here, [S*] is concentration for all the protein sub- 
stances in the phlegm agglomerate. Compared with pro- 
tein concentration [So] of the SARS coronavirus existing 
alone, there satisfies such a relation, i.e. [S*] = «[So], 
where « is ratio for the volume of all protein substances 
to that of the SARS coronavirus. Since the virus was 
contained in the phlegm or saliva, if only 90% of the 
protein substance become denatured, can it be said the 
virus was killed. Assuming the SARS coronavirus only 
takes up 10% volume of the total phlegm or saliva, one 
has « = 10. Then the increased survival time for the virus 
will be In 10 times of its original value. For example, an 
original survival time of 5 h will now increase to 
5ln10 = 11.5 h. Therefore, cleaning away the environ- 
mentally existing phlegm or saliva is very necessary for 
the successful control of the SARS infection. 


6.7 

Implementation of rate equation theory to therapy 

and detection of SARS 

The above study demonstrates for the first time that, the 
rate equation has flexible adaptability to characterize the 
survival time of the SARS coronavirus. Once the activation 
energy and frequency factor in a specific environment are 
available, one can perform predictions on the survival 
behavior of the SARS coronavirus in a wide variety of 
surrounding situations. The present theory appears useful 
not only for better understanding the spreading and sur- 
vival mechanisms of the SARS coronavirus, it also suggests 
an important clue to develop efficient drugs, bacterin and 
clinical therapy. For example, if certain strategy capable of 
varying the activation energy, frequency factor and tem- 
perature of the SARS coronavirus can be found, it is very 
possible to set up a useful way for fighting the SARS. In 
fact, such effort has been reflected in the current SARS 
research activities. Compared to the ways of developing 
bacterin and drug, physical therapy appears much direct 
in its output results. For example, by testing the killing 
effect of a group of various physical factors to the in vitro 
and in vivo SARS coronavirus, it is possible to quickly 
screen out useful engineering approaches for treating the 
SARS. The currently available approaches worthy of trying 
are as follows [31]: Based on the temperature sensitivity of 
the SARS coronavirus, one can develop various heating 
methods for the hyperthermia and disinfection purposes; 
While from the non-thermal effects produced by the 
electromagnetic wave (static or dynamic electric stimula- 
tion, radio frequency), the optical wave (laser or violet), 
acid and alkali environment (by varying the pH value), the 
mass diffusion (high concentration oxygen or human 
tolerable disinfector), or delivered target substance (drug 
or bacterin) etc., one can also find useful ways to intervene 
the rate process of the SARS coronavirus and then suc- 
cessfully prevent or cure it. Similarly, if certain specific 
substance was found to react with the samples containing 
SARS coronavirus and produce an easily detectable 
physical chemical mass, one can establish a quick diag- 
nostic method on the SARS in clinics. 


7 

Conclusions 

It is a good way to quantitatively characterize the 
denaturation or replication behavior of SARS coronavi- 
rus through the rate process theory. Based on the limited 
experimental data on survival time, the activation energy 
and frequency factor for the SARS coronavirus in a 
specific environment can be estimated and then applied 
to predict more information inaccessible through expe- 
riment. The currently available experimental data for the 
survival times of SARS coronavirus is very limited, and 
the measurements only reflect certain specific cases. It is 
very necessary to perform systematic tests on the sur- 
vival times of SARS coronavirus subject to electric, 
magnetic, acoustic, optical, thermal or acid-alkali effects 
with different strengths. Only in this way, can the data 
bank for the activation and factor be established. The 
current data for the survival times of the SARS coro- 
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navirus provided by different labs appears a little dis- 
perse from each other. Their accuracy is in need of 
re-evaluation. One possible way can be that, the same SARS 
coronavirus from the only source should be distributed to 
multiple labs over the world for independent tests on the 
survival time of the sample under the same experimental 
conditions strictly defined. In this way, deviation from each 
test can be tested and thus allows for a comprehensive 
understanding on the SARS coronavirus. Such strategy will 
also be helpful for testing the variability of SARS corona- 
virus. The causes to kill the SARS coronavirus can be at- 
tributed to the effects of changing the activation energy and 
frequency factor. Any drug or engineered way capable of 
inducing a death of SARS coronavirus will be beneficial for 
the treatment of SARS. Therefore, it is critical to perform 
experimental study on the survival or replication of SARS 
coronavirus either in vivo or vitro. Through a compre- 
hensive analysis on the effects of various physical and 
chemical factors, one can establish the possible drug or 
engineering way. Although the present study is mainly for 
the SARS coronavirus, behaviors for any general virus to 
attack its host cell are approximately similar to each other 
and therefore can be described by the rate process theory. 
Analysis given in this paper will be beneficial in quantify- 
ing the physical or chemical behavior of SARS coronavirus. 
Up to now, little is known in this area. The present study 
provides a new theoretical foundation. 
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